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Abstract 

The polarization of cosmic microwave background (CMB) can be 
used to constrain cosmological birefringence, the rotation of the lin¬ 
ear polarization of CMB photons potentially induced by parity vi¬ 
olating physics beyond the standard model. This effect produces 
non-null CMB cross correlations between temperature and B mode- 
polarization, and between E- and B-mode polarization. Both cross¬ 
correlations are otherwise null in the standard cosmological model. 

We use the recently released 2015 Planck likelihood in combina¬ 
tion with the Bicep2/Keck/Planck (BKP) likelihood to constrain the 
birefringence angle a. Our findings, that are compatible with no de¬ 
tection, read a = 0.0° ± 1.3° (stat) ± 1° (sys) for Planck data and 
a = 0.30° ± 0.27° (stat) ± l°(sys) for BKP data. 

We finally forecast the expected improvements over present con¬ 
straints when the Planck BB, TB and EB spectra at high I will be 
included in the analysis. 

1 Introduction 

It has become customary to use polarized cosmic microwave back¬ 
ground (CMB) data to constrain the cosmic birefringence effect [H [2l 
|3], i.e. the in vacuo rotation of the photon polarization direction dur¬ 
ing propagation [3]. In general, such effect, which is naturally param¬ 
eterized by an angle a, results in a mixing between Q and U Stokes 
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parameters that produces non-null CMB cross correlations between 
temperature and B mode-polarization, and between E- and B-mode 
polarization. Since these correlations are expected to be null under 
the parity conserving assumptions that is beneath the standard cos¬ 
mological model, cosmic birefringence is a tracer of parity violating 
physics. 

Several theoretical models exhibit cosmological birefringence. Any 
coupling of the electromagnetic (EM) field with axion-like particles 
[5] or a quintessence field [6] is expected to induce the effect, which 
can be also driven by quantum-gravity terms [7] or Chern-Simons type 
interactions [3] in the EM Lagrangian. Each particular model predicts 
a specific dependence of the CMB angular power spectra (henceforth 
APS) on a, which can be complicated. In this paper, we restrict to the 
simplest case of constant a, for which the effect can be parametrized 
as [Hl|9] 


^TE^obs 

= Cj^ cos(2a), 

(1) 

^TB,obs 

II 

do 

(2) 

^EE,obs 

= Cf^ cos^(2a) -|- sin^(2a), 

( 3 ) 

^BB,obs 

= cos^(2a) -|- sin^(2a), 

( 4 ) 

^EB^obs 

= A - Cf^) sin(4a), 

( 5 ) 


with and being the observed and the primordial (i.e. 

unrotated) APS for the XY spectrum {X, Y = T, E or B), i.e. the 
one that would arise in absence of birefringence. Note that we set 
the primordial TB and EB spectra to zero for the sake of simplic¬ 
ity, thereby assuming parity violation effect plays a negligible role up 
to CMB photon decoupling (this choice excludes e.g. chiral gravity 
theories). 

The existing limits on a, as obtained from CMB experiments [H 
m El uni El [12], are shown in Fig. (see first seven estimates). In 
short, almost all constraints are well compatible with a null effect, 
and the best available statistical errors are of order of 0.2° with a 
systematic uncertainty estimated at the level of 0.5°. For comparison, 
still in Fig. we anticipate the constraints we obtain from Planck 
and Bicep2/Keck/Planck (henceforth BKP) data in this paper (see 
the last two constraints). As discussed below the Planck constraints 
do not include B information at i > 29 while this is included in BKP 
combination. 

This paper is organized as follows: in Section we provide a de¬ 
scription of the performed analysis; in Section[^we report the datasets 
employed and the associated likelihood functions; in Section we 
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Figure 1: Constraints on a from several CMB experiments HI El El DEI El [12]. 
Thin black error bars are for statistical uncertainties at 68% C.L.. Thick blue 
error bars are obtained by summing linearly the statistical and systematic 
uncertainties. The statistical uncertainty of BOOM03 do already contain a 
contribution from systematic effect. The last two constraints are derived in 
this paper. 


3 
















present the main constraints on a. sampling on all standard ACDM 
parameters and as well on the tensor-to-scalar ratio, r. We also discuss 
the role of systematic uncertainties. We finally draw our conclusions 
in Section providing also forecast for the statistical sensitivity that 
Planck data can reach when BB, TB and EB spectra will be included 
at high <1 in the Planck likelihood. 


2 Description of the analysis 

CMB polarization arises at two distinct cosmological times: the re¬ 
combination epoch (z ~ 1100) and the reionization era at 2 ; ~ 10 
m- When the CMB field is expanded in spherical harmonics, the 
first signal mostly shows up at high multipoles, since polarization is 
generated through a causal process and the Hubble horizon at last 
scattering only subtends a degree sized angle. The later reionization 
of the cosmic fluid at lower redshift impacts instead the lowest mul- 
tipoles. These two regimes need to be properly treated when probing 
for cosmological birefringence. It is also interesting to separate the 
constraints on a between high and low i as these can be ascribed to 
different epochs and, hence, physical conditions. Therefore, following 
|14j . we perform a Monte Carlo Markov chain (MCMC) analysis mak¬ 
ing use of the publicly available code cosmomc [mnB] to constrain a 
by considering two multipole range^ 2 < i < 29 and 30 < £ < 2500. 

In practice, we first sample at low £ on a, on the optical depth r 
and on the scalar amplitude Mg, while fixing the remaining ACDM 
cosmological parameters to the Planck best-fit model derived from 
temperature and polarization (TT, TE, EE). We then focus on the 
high £, fixing r to the best-fit value obtained in the former analysis 
and let the remaining parameters to vary. We call the above combina¬ 
tion case A. We report separately the constraints on a obtained by a 
single run on all multipoles. This is called case B in what follows and 
it is provided for comparison purposes. Moreover, we consider a third 
case, C, where we sample at the same time on two separate angles, 
namely aiow and a^igh, which rotate the APS independently in the 
two multipole regimes. 

These three cases above are performed considering only scalar per¬ 
turbations (i.e. ACDM-l-a) and also considering tensor perturbations 
(i.e. ACDM-|-r-|-Q;). See Tablefor a summary of all cases considered. 

^The transition multipole between low and high i is taken to be 29 following the low 
and high transition of the Planck likelihood (see also Section]^. The fine choice of this 
multipole does not impact significantly the results as long as it occurs in a region where 
the polarization signal is noise dominated m- 
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Table 1: Summary of the considered cases. See text. 


case 

description 

color code 

A 

First run at low £ on a, Ag and r 
(other parameters hxed to Planck best £t) 

light blue 


Second run at high £ on ACDM(-|-r)-|-a 
with T fixed to the first run best fit value 

blue 

B 

Single run on ACDM(-|-r)-|-Q; in all the £ 
range 

red 

C 

Single run on ACDM(-|-r)-|-Q;jo^-fQ:/ijg/i 
where aiow (dashed) and ahigh (solid) are 
defined in the low and high £ range 
respectively 

green 


Gravitational leasing induces non zero TB and EB correlations 
through E to B mixing, an effect potentially observable at the sensi¬ 
tivity level of the Planck and BKP dataset. To disentangle it from 
cosmological birefringence, one should in principle model the latter 
within the CMB time evolution from last scattering to us. However, 
it can be shown [18] that for the model given in Eqs. ([T]j^, one can 
safely separate the birefringence and leasing contributions. 


3 Datasets and likelihoods 

We have used the Planck and BKP publicly available likelihood^ 

Here we list the several terms entering the likelihood functions [I71II9|. 

• Planck lowl. Planck low i temperature likelihood. Based on 
the Planck 2015 low resolution CMB anisotropy map (Com¬ 
mander, see [20jh 

• Planck lowTEB. Planck low t temperature and polarization 
likelihood. Based on temperature (Commander) and polariza¬ 
tion (70GHz) maps. 

• Planck TT,TE,EE. Planck high t' temperature and polarization 
Likelihood. Based on TT, TE and EE cross-spectra. 

• BKP BB. Bicep2/Keck likelihood with 353GHz Planck channel. 

Based on the hrst hve band powers of the BB spectra m- 

^For Planck see http://pla.esac.esa.int/pla/ and for BKP see 
http: //bicepkeck.org/bkp _2015_release.html 
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Figure 2: Left panel; Constraints on a from the three considered cases sum¬ 
marized in Table [1} Blue curves are for case A, red curves for case B and 
green curves for case C (solid green for £ G [30, 2500] and dashed green for 
£ G [2,29]). Right panel: zoom of the left panel. 



• BKP EE,BB. Bicep2/Keck likelihood with 353GHz Planck chan¬ 
nel. Based on the first five band powers of the BB, EE and EB 
spectra. 


From Eqs. (I][5) it is clear that we are sensitive to the sign of a 
only when we use TB and EB information. This happens only when 
we employ Planck lowTEBandBKP EE,BB. Therefore, when the latter 
are not considered, in fact our constraints on a will be on its absolute 
value. 

For each considered case, we employ the following likelihoods: 


• Case A. First run Planck lowTEB. Second run Planck lowl -|- 
Planck TT,TE,EE. 

• Case B. Single run with Planck lowTEB + Planck TT,TE,EE. 

• Case C. Single run with Planck lowTEB -|- Planck TT,TE,EE. 


Ah the above cases are analyzed for both ACDM-Pa and ACDM-|-r-|-a. 
Moreover, in the latter case we consider the likelihoods BKP BB or BKP 
EE,BB either alone (with minimum multipole £ ~ 50) or in combination 
with Planck. 


4 Results 
4.1 ACDM+ct 

The ACDM-|-a model is analyzed considering Planck data alone. 
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Figure 3: Posterior distributions of a, r and Ag for the low i run of case A. 
Blue curves consider temperature and polarization at low ^ while light blue 
take into account polarization at low £ and temperature at all 

The posterior distributions for a are shown in Fig. and the corre¬ 
sponding constraints are reported in Table The latter only account 
for statistical uncertainty (68% C.L.). The companion “low-£ run” of 
case A is displayed in Fig. i Fig. i shows that all the standard 
ACDM cosmological parameters are very stable against the inclusion 
of a in the MCMC and Fig. shows how the same parameters are 
very stable in all three considered cases. 

All the Planck constraints are compatible with zero even con¬ 
sidering only the statistical uncertainty. Still without accounting for 
systematic uncertainty, the best constraint for Planck alone is a = 

0.0° ± 1.3°. The standard error will improve significantly in a future 
version of the high ^ likelihood which will employ the TB and EB 
spectra. Forecasts are provided below. 

4.2 ACDM+r+ct 

The ACDM-|-r-|-Q; model is analyzed with the following data set com¬ 
binations: Planck alone, BKP and Planck plus BKP. Since all data 
in the BKP dataset come from the small {fsky — 1%) region observed 
by Bicep, we neglect the bias caused by counting twice the Planck 
dataset in this region, for the latter combination. As for the ACDM-Pa 
model, all the estimates are reported in Table 

• Planck data alone. The posterior distribution functions of a 
are displayed in Fig. In Fig. we show the low i compan¬ 
ion run of case A. The posterior distribution functions of the six 
standard cosmological parameters and r are plotted in Fig. 

All the cases are very stable and consistent. The inclusion of 
OL does not produce any shift in the ACDM-|-r parameters. As 
for the ACDM-|-a model, even in the ACDM-|-r-|-Q; model all the 



a\owi [°] T ln(10^°As) 
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Figure 4: Comparison of the base ACDM model parameters and the bire¬ 
fringence angle a from the analysis over the 2 < i < 2500 range (case B). All 
the 2D contour plots are very stable against the inclusion of a in the MCMC 
sampling. We conclude that a is largely uncorrelated with the other ACDM 
parameters. 






















Table 2: Constraints of a for the posterior distributions shown in Fig. 
Units: deg. P stands for PLANCK, BKP EE and BKP EE BB for Bi- 
cep2/Keck/Planck BB and EE,BB,EB data, unc = unconstrained, a cor¬ 
responds to the case in which the six ACDM parameters are hxed to the 
Planck best £t. 


case 

^low 

C^high 

a 

model 

ACDM 

data 

A 

c: + 12 

^-15 

0.0 ± 1.4 

- 

+a 

P 

B 

- 

- 

-0.03 ± 1.35 

+a 

P 

C 


0.0 ± 1.3 

- 


P 

A 


0.0 ± 1.4 

- 

±a±r 

P 

B 

- 

- 

0.00 ± 1.33 

±a±r 

P 

C 

0.6t?3% 

0.0 ± 1.3 

- 

±a±r 

P 

A 


0.0 ± 1.4 

- 

±a±r 

P±BKP BB 

B 

- 

- 

0.02tjij 

±a±r 

P ± BKP BB 

C 

0.5^?o®o 

0.0 ± 1.4 

- 

±a±r 

P ±BKP BB 

A 

4tl"5 

0.32 ±0.27 

- 

±a±r 

P±BKP BB EE 

B 

- 

- 

0.32 ±0.26 

±a±r 

P ± BKP BB EE 

C 

n q+9-9 

'^•^-13.6 

0.32 ±0.26 

- 

±a±r 

P ±BKP BB EE 

C* 

unc 

0.0 ± 1.5 

- 

±a±r 

BKP BB 

C* 

unc 

0.30 ±0.27 

- 

±a±r 

BKP BB EE 


- ACDM,Planck TT,TE,EE+lowTEB - ACDM+ah«, r from lowTEB,Planck TT,TE,EE+lowl 

- ACDM+a,Planck TT,TE,EE+lowTEB - ACDM+a/o^ + a^/r,Planck TT,TE,EE+lowTEB 



Figure 5: Marginalized posterior distributions for the six main cosmological 
parameters for the ACDM when jointly sampled with a for the three consid¬ 
ered cases. All the parameters are very compatible among the three cases. 
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Figure 6: Left panel: Constraints on a from the three considered cases taken 
into account in ACDM+r+a. Blue curves are for case A, red curves for case 
B and green curves for case C (solid green for 30 — 2500 and dashed green for 
2 — 29). Right panel is just a zoom of the right panel. Only Planck data 
are employed. 
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Figure 7: Low i run of case A. In ACDM+r+a. 
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ACDM+r,Planck TT,TE,EE+lowTEB 
ACDM + r + a.Planck TT.TE.EE+lowTEB 


ACDM+r + a^ii, t from lowTEB,Planck TT,TE,EE+lowl 
ACDM + r + aioe + Planck TT,TE,EE+lowTEB 



Figure 8: Marginalized posterior distributions for the seven main cosmologi¬ 
cal parameters for the ACDM-l-r (blue) and ACDM-|-r-f-a (red) models from 
the analysis over the 2 < i < 2500 range. It is clear that the parameters are 
very stable against the inclusion of a. 

Planck constraints on a are compatible with zero and numeri¬ 
cally very similar. 

• Planck and BKP data. The probability distribution functions 
of a are shown in Fig. Blue, red and green solid (dotted- 
dashed) curves are for Planck -|- BKP BB (BKP EE BB) data 
in the cases A, B and C respectively. Cyan curves are for BKP 
data alone (solid for BKP BB and dotted-dashed for BKP EE 
BB). The corresponding constraints are reported in Table 
When Planck data are employed in combination with BKP BB 
we notice a change in the shape of the probability distribution 
function and a bimodal distribution appears. This is due to the 
fact that the sign of ahigh is not constrained when EB and TB 
spectra are neglected. However the compatibility with zero is 
at the same level as in the Planck data alone analysis. When 
BKP EE BB data are taken into account in combination with 
Planck, due to the presence of the EB spectrum, the bimodal 
pattern disappears and a better constraint of a is obtained, i.e. 
a = 0.32° ± 0.26° (again, only statistical). In this case the BKP 
data dominates the combination and the impact of Planck data 
is subleading since what counts more is the EB spectrum. 
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Figure 9: Probability distribution function for a. Blue, red and green solid 
(dotted-dashed) curves are for Planck -|- BKP BB (BKP EE BB) data in 
the cases A, B and C respectively. Cyan curves are for BKP data alone (solid 
for BKP BB and dotted-dashed for BKP EE BB). 
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4.3 Uncertainties from systematic effects 

The most relevant untreated systematic effect that impacts birefrin¬ 
gence constraints is a miscalibration of the detectors’ optimal polar¬ 
ization angle, i.e. the angle, measured in the telescope reference frame, 
at which a detector is maximally sensitive to a vertical linear polariza¬ 
tion m- In the case of Planck, ground based measurements cannot 
be assumed to be representative of in-flight conditions at face value, 
but rather need to be extrapolated since thermoelastic effects arising 
in operations can play a role. In-flight measurements are complicated 
by the scarsity of linearly polarized sources that are bright enough, 
though the Crab nebula has been attempted as a calibrator m- For 
Planck, we use the estimate of the error in the polarization orien¬ 
tation evaluated in [23j, where it is given at 1°. The original Bicep2 
analysis [22] discusses how the TB and EB spectra are rotated to min¬ 
imize E to B mode mixing, by imposing that their residuals from zero 
are minimized. This process destroys of course any detection power for 
birefringence below this final rotation angle, which is estimated at 1° 
and can thus be taken as a measure of the final systematic budget. In 
the BKP analysis, the treatment of the Bicep and Keck power spectra 
is identical to that of Bicep2. We therefore assume a 1° systematic un¬ 
certainty of the instrumental polarization angle for the BKP dataset 
as wel0 The same assumption applies to the joint analysis of BKP 
and Planck. 


5 Conclusions 

We have constrained the birefringence angle a employing Planck 
and BKP data either alone or in combination. All our results are well 
compatible with no detection even without accounting for systematic 
uncertainty. We find a = 0.0° ± 1.3° (stat) ± 1° (sys) for Planck 
data alone. Similar results are found when we also consider the BB 
spectrum of BKP data. The inclusion of the EB spectrum in the BKP 
likelihood makes the constraint on a much tighter, namely ol = 0.32° ± 

0.26° (stat.) ± l°(sys.). In the latter case the inclusion of Planck data 
brings little constraining power. This is expected since the current 
version of the Planck Likelihood does not contain TB and EB spectra 
at high i. It is possible to forecast that including BB, TB and EB 
spectra in the high i Planck likelihood and considering the range of 
multipoles, i.e. [30-2500], the level of the statistical uncertainty on a 
can be as low as ~ 0.03°. This is obtained considering the nominal 

^Notice that for Keck the rotation angle is smaller, i.e. about 0.5° degrees [53]. We 
assume therefore that the error is dominated by Bicep. 
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sensitivity of the 143 GHz channel and an exact likelihood approach, 
see e.g. [2H]. The obtained value is compatible with what forecasted in 
the Planck bluebook [26] considering the longer integration time and 
the improved sensitivity of the nsed channel. Of course, the quest for 
reducing the statistical errors would eventually prove meaningless if it 
is not accompanied by a parallel effort to understand and reduce the 
systematic error budget affecting polarization rotation constraints. 
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